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Abstract 

The effect of orbital magnetism on the chemical bonding of lateral, two-dimensional 
artificial molecules is studied in the case of a 2e double quantum dot (artificial molec- 
ular hydrogen). It is found that a perpendicular magnetic field reduces the coupling 
(tunneling) between the individual dots and, for sufficiently high values, it leads to 
complete dissociation of the artificial molecule. The method used is building on 
Lowdin's work on Projection Operators in Quantum Chemistry; it is a spin-and- 
space unrestricted Hartree-Fock method in conjunction with the companion step of 
the restoration of spin and space symmetries via Projection Techniques (when such 
symmetries are broken). This method is able to describe the full range of couplings 
in two-dimensional double quantum dots, from the strong-coupling regime exhibiting 
delocalized molecular orbitals to the weak-coupling and dissociation regimes associ- 
ated with a Generalized Valence Bond combination of atomic-type orbitals localized 
on the individual dots. 
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1 Introduction 



Two-dimensional (2D) Quantum Dots (QD's) [[T], U are usually referred to as artificial 
atoms, a term suggestive of strong similarities between these manmade nanodevices and 
the physical behavior of natural atoms. As a result, in the last few years, an intensive the- 
oretical effort [2—13] has been devoted towards the elucidation of the appropriate analogies 
and/or differences. Recently, we showed [10—12] that in the absence of a magnetic field the 
most promising analogies are mainly found outside the confines of the central-field approxi- 
mation underlying the Independent-Particle Model (IPM) and the ensuing physical picture 
of electronic shells and the Aufbau Principle. Indeed, as a result of the lower electronic den- 
sities in QD's, strong e — e correlations can lead (as a function of the ratio Rw between the 
interelectron repulsion and the zero-point kinetic energy) to a drastically different physical 
regime, where the electrons become localized, arranging themselves in concentric geometric 
shells and forming electron molecules (referred to also as Wigner molecules in analogy to 
Wigner crystallization |14| in infinite media). In this context, it was found |ll], that 
the proper analogy for the particular case of a 2e QD is the collective-motion picture rem- 
iniscent of the fleeting and rather exotic phenomena of the doubly-excited natural helium 
atom, where the emergence of a "floppy" trimeric molecule (consisting of the two localized 
electrons and the heavy a-particle nucleus) has been well established [p~5] , [16]. 

A natural extension [10,13,17—23] of this theoretical effort has also developed in the 
direction of 2D QD Molecules (QDM's, often referred to as artificial molecules), aiming 
at elucidating the analogies and differences between such artificially fabricated |24|, |25 



molecular nanostrustures and the natural molecules. (Depending on the arrangement of 
the individual dots, two classes of QDM's can be distinguished: lateral [10,13,17—21,24] 
and vertical |22|, |23|, |25[ ones.) 



In a previous paper [21] and for the case of zero magnetic field (field-free case), we 
addressed the interplay of coupling and dissociation in lateral QDM's. We showed that this 
interplay relates directly to the nature of the coupling in the artificial molecules, and in 
particular to the question whether such coupling can be described by the Molecular Orbital 
(MO) Theory or the Valence Bond (VB) Theory in analogy with the chemical bond in 
natural molecules. 

A major attraction of QD's and QDM's is the fact that, due to their larger size, orbital 
magnetic effects become important for magnetic-field values easily attainable in the labo- 
ratory. This contrasts with the case of natural atoms and molecules for which magnetic 
fields of sufficient strength for the production of phenomena related to orbital magnetism 
are known to occur only in astrophysical environments, e.g., on the surfaces of neutron stars 



26], [27[] . In this paper, we study the effect of orbital magnetism on the interdot coupling in 
a lateral, two-electron double quantum dot. This nanodevice represents an artificial analog 
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to the natural hydrogen molecule and can be denoted as H 2 -QDM. It has been suggested 
|H| that it can function as the elemental two-qubit logic gate in quantum computing. 

We find that the interplay of the MO versus the VB description provides the proper 
framework for understanding the influence of orbital magnetism on the chemical bonding 
of the H2-QDM. In particular, we show that a perpendicular magnetic field reduces the 
coupling between the individual dots and, for sufficiently high values, it leads to the dis- 
sociation of the artificial molecule. As a result, in addition to the obvious parameters of 
interdot barrier height and interdot separation, the magnetic field supplies a third variable 
able to induce dissociation and thus to control the strength of the interdot coupling. 

Our approach is twofold. As a first step, we utilize a self-consistent-field theory which 
can go beyond the MO approximation, namely the spin-and-space unrestricted Hartree- 
Fock (sS-UHF), which was introduced by us |TD|, |TT| for the description of the many-body 
problem of both single [[K], and molecular JTOj] QD's. The equations used are given in 
Ref. Iffijjl , where they are simply referred to as unrestricted Hartree-Fock (UHF); the addi- 
tional sS labeling employed by us emphasizes the range of possible symmetry unrestrictions 
in the solutions of these equations. In particular, the sS-UHF differs from the more fa- 
miliar restricted HF (RHF) in two ways: (i) it relaxes the double- occupancy requirement — 
namely, it employs different spatial orbitals for the two different (i.e., the up and down) 
spin directions [thus the designation "spin (s) unresricted"], and (ii) it relaxes the require- 
ment that the electron orbitals be constrained by the symmetry of the external confining field 
[thus the designation "space (S) unrestricted"]. Since it is a general property of the HF 
equations to preserve at each iteration step the symmetries of the many-body hamiltonian 
(whenever they happen to be present in the HF electron density), the input trial density 
at the initial step must be constructed in such a way as to a priori reflect the relaxation 
of the two requirements mentioned above. Observe further that, in order to describe elec- 
tron localization, the sS-UHF employs fully the fact that all iV (where N is the number of 
electrons) orbital-dependent effective (mean-field) HF potentials can be different from each 
other. 

We remark that within the terminology adopted here, the simple designation Hartree- 
Fock (HF) in the literature most often refers to our restricted HF (RHF), in particular 
in atomic physics and the physics of the homogeneous electron gas. In nuclear physics, 
however, the simple designation HF most often refers to a space (S)-UHF. The simply 
designated unrestricted Hartree-Fock (UHF) as used in Chemistry (e.g., in calculations of 
open shell molecules) corresponds most often to our s-UHF (but not simultaneously space 
unrestricted HF). 

As a second step, we will show that, in conjunction with Lowdin's spin-projection tech- 
nique p0| , the solutions with broken space symmetry allowed in QDM's by the sS-UHF 



provide a natural vehicle for formulating a Generalized Valence Bond (GVB) theory (see 
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below) able to describe the chemical bonding in artificial molecules in both the case of an 
applied magnetic field and the field-free case. 



2 The two-center-oscillator confining potential 

In the 2D two-center-oscillatorQ (TCO), the single-particle levels associated with the con- 
fining potential of the artificial molecule are determined by the single-particle hamiltonian 



32, 33 
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tt rri | - L *22,- L *2/2 

H = T + -m u xk x + -m u yk y k 

+ V neck (y) + h k + • s , (1) 

where y' k = y—y k with k — 1 for y < (left) and k = 2 for y > (right), and the h k s control 
the relative well-depth, thus allowing studies of hetero-QDM's. x denotes the coordinate 
perpendicular to the interdot axis (y). T = (p — eA/c) 2 /2m*, with A = 0.5(— By, Bx, 0), 
and the last term in Eq. (|I|) is the Zeeman interaction with g* being the effective g factor, 
fis the Bohr magneton, and s the spin of an individual electron. Here we limit ourselves 
to systems with h\o x \ = hu X 2 = htv x . The most general shapes described by H are two 
semiellipses connected by a smooth neck [V neck (y)]. yi < and y 2 > are the centers of 
these semiellipses, d = yi~y\ is the interdot distance, and m* is the effective electron mass. 

For the smooth neck, we use V neck {y) = lm*Uy k [c k y' k 3 + d k y k }9(\y\ - \y k \), where d(u) = 
for u > and 9{u) = 1 for u < 0. The four constants c k and d k can be expressed via two 
parameters, as follows: (— l) k c k = (2 — 4e b k )/y k and d k — (1 — 3e b k )/yl, where the barrier- 
control parameters e b k = (VJ, — h k )/V ok are related to the actual (controlable) height of the 
bare barrier (V&) between the two QD's, and Vo k = m*ujy k y k /2 (for hi = h 2 , Voi = V02 = Vo). 

The single-particle levels of H, including an external perpendicular magnetic field B, 
are obtained by numerical diagonalization in a (variable-with-separation) basis consisting 
of the eigenstates of the auxiliary hamiltonian: 

d 2 1 1 

tt P l - L *22|- L *2/2,7, /r>\ 

= 2^ 2 m UxX 2 m UykVk k ' ( ' 

This eigenvalue problem is separable in x and y, i.e., the wave functions are written as 
&mu(%, y) = X m (x)Y v (y). The solutions for X m (x) are those of a one-dimensional os- 
cillator, and for Y u (y) they can be expressed through the parabolic cylinder functions 



32], |33| U[a k , (-l) fe ^ fc ], where £ k = y' k ^/2m*uj yk /h, a k = (-E y + h k )/(fkv yk ), and E : 



u 



1 A 3D magnetic-field- free version of the TCO has been used in the description of fission in metal clusters 
52| and atomic nuclei IB3|]. 
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{y + 0.5) fitd y i + hi denotes the y-eigenvalues. The matching conditions at y = for the left 
and right domains yield the y-eigenvalues and the eigenf unctions Y u (y) (m is integer and v 
is in general real). 

In this paper, we will limit ourselves to symmetric (homopolar) QDM's, i.e., hw x = 
fojjyi = fojjyi = hujQ, with equal well-depths of the left and right dots, i.e., hi = h 2 = 0. In 
all cases, we will use hwo = 5 meV and m* = 0.067m e (this effective-mass value corresponds 
to GaAs). 



3 The Many-Body Hamiltonian 

The many-body hamiltonian 7i for a dimeric QDM comprising N electrons can be expressed 
as a sum of the single-particle part H{i) defined in Eq. ([I]) and the two-particle interelectron 
Coulomb repulsion, 

N N N 2 

i=i i=i j>i j 

where k is the dielectric constant and r^- denotes the relative distance between the i and j 
electrons. 

As we mentioned in the introduction, we will use the sS-UHF method for determining at 
a first level an approximate solution of the many-body problem specified by the hamiltonian 
(0). The sS-UHF equations are solved in the Pople-Nesbet-Roothaan formalism [28| using 



the interdot-distance adjustable basis formed with the eigenfunctions & mu (x, y) of the TCO 
defined in section 2. 

As we will explicitly illustrate in section 4.2 for the case of the H2-QDM, the next step 
in improving the sS-UHF solution involves the use of Projection Techniques in relation to 
the UHF single Slater determinant. 



4 Artificial molecular hydrogen (H2-QDM) in a mag- 
netic field: A Generalized Valence Bond approach 

4.1 The sS-UHF description 

As an introductory example to the process of symmetry breaking in HF, we consider in this 
subsection the field-free (B = 0) case of H2-QDM with k = 20 (this value is an intermediate 
one to the three different values of k that will be considered below in the case of an applied 
magnetic field). Fig. 1 displays the RHF and sS-UHF results for the P = N a — N@ = 
case (singlet) and for an interdot distance d = 30 nm and an interdot barrier V& = 4.95 
meV (a and j3 denote up and down spins, respectively). In the RHF (Fig. 1, left), both 
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the spin-up and spin-down electrons occupy the same bonding (a g ) molecular orbital. In 
contrast, the sS-UHF results exhibit breaking of the spatial reflection symmetry; namely, 
the spin-up electron occupies an optimized?] Is atomic-like orbital (AO) in the left QD, while 
the spin down electron occupies the corresponding Is' AO in the right QD. Concerning the 
total energies, the RHF yields Erhf(P = 0) = 13.68 meV, while the sS-UHF energy is 
E s suhf(P = 0) = 12.83 representing a gain in energy of 0.85 meV. Since the energy of the 
triplet is E sUHF (P = 2) = E sS uhf(P — 2) = 13.01 meV, the sS-UHF singlet conforms to 
the requirement |34] that for two electrons at zero magnetic field the singlet is always the 
ground state; on the other hand the RHF MO solution fails in this respect. 



4.2 Projected wave function and restoration of the broken sym- 
metry 

To make further progress, we utilize the spin projection technique to restore the broken 
symmetry of the sS-UHF determinant (henceforth we will drop the prefix sS when referring 
to the sS-UHF determinant), 



V2^ UHF (1,2) = 



«( ri )a(l) v( ri )P(l 
u(r 2 )a(2) v{v 2 )(3{2) 



= \u(l)v(2) > , (4) 



where u(r) and v(r) are the Is (left) and Is' (right) localized orbitals of the sS-UHF solution. 
An example of such orbitals for the field-free case are displayed in the right column of Fig. 
1. Similar localized orbitals appear also in the B^O case, so that in general the functions 
u(r) and v (r) are complex, a and (3 denote the up and down spin functions, respectively. 
In Eq. (f|) we also define a compact notation for the ^uhf determinant, where a bar over 
a space orbital denotes a spin-down electron; absence of a bar denotes a spin-up electron. 

^uhf{1, 2) is an eigenstate of the projection S z of the total spin S = Si + s 2 , but not 
of S 2 . One can generate a many-body wave function which is an eigenstate of S 2 with 
eigenvalue s(s + l)h 2 by applying the following projection operator introduced by Lowdin 



30, 31 



P = n 82 - + l)h2 (5) 



where the index s' runs over the quantum numbers of S . 



2 The optimized orbitals are anisotropic (i.e., non-circular ly symmetric) reflecting polarization effects 
due to the electronic interdot interaction. 
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FIGURE 1 



Figure 1: Lateral H 2 -QDM at zero magnetic field: Occupied orbitals (modulus square, 
bottom half) and total charge (CD) and spin (SD) densities (top half) for the P = spin 
unpolarized case. Left column: RHF results. Right column: sS-UHF results exhibiting 
a breaking of the space symmetry. The numbers displayed with each orbital are their 
eigenenergies in meV, while the up and down arrows indicate an electron with an up or 
down spin. The numbers displayed with the charge densities are the total energies in 
meV. Unlike the RHF case, the spin density of the sS-UHF exhibits a well developed spin 
density wave. Distances are in nm and the electron densities in 10~ 4 nm~ 2 . The choice of 
parameters is: m* = 0.067m e , Hcuq = 5 meV, d = 30 nm, V& = 4.95 meV, k = 20. 
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The result of S 2 on any UHF determinant can be calculated with the help of the ex- 
pression, 



s 2 <$> UHF = h 2 



(N a -Np) 2 /4 + N/2 + J2™ 



i<j 



&UHF , (6) 



where the operator tUy interchanges the spins of electrons % and j provided that their 
spins are different; N a and N@ denote the number of spin-up and spin-down electrons, 
respectively, while N denotes the total number of electrons. 

For the singlet magnetic state of two electrons (N = 2), one has N a = Np = 1, S z = 0, 
and S 2 has only the two quantum numbers s = and s = 1. As a result, 

uhf{^, 2) 

= \u(l)v(2) > - \u{l)v{2) > . (7) 

In contrast to the single-determinantal wave functions of the RHF and sS-UHF methods, the 
projected many-body wave function (|7|) is a linear superposition of two Slater determinants, 
and thus represents a corrective step beyond the mean-field approximation. 
Expanding the determinants in Eq. (|7|), one finds the equivalent expression 

2P * VHF (1, 2) = (u( ri Mr 2 ) + u (r 2 )«( ri ))x(0, 0) , (8) 

where the spin eigenfunction is given by 

X (s = 0;S Z = 0) = (a(l)/3(2) - a(2)/3(l))/ v / 2 . (9) 

Eq. (§) has the form of a Heitler-London (HL) or valence bondP] |36], ^7] wave function 



for the singlet magnetic state. However, unlike the HL scheme which uses the orbitals 
4>l(t) and <Ar(i") of the separated (left and right) atoms,f] expression (JD employs the sS- 

3 The early empirical electronic model of valence was primarily developed by G.N. Lewis who introduced 
a symbolism where an electron was represented by a dot, e.g., H:H, with a dot between the atomic symbols 
denoting a shared electron. Later in 1927 Heitler and London formulated the first quantum mechanical 
theory of the pair-electron bond for the case of the hydrogen molecule. The theory was subsequently 
developed by Pauling and others in the 1930's into the modern theory of the chemical bond called the 
Valence Bond Theory. 



4 Refs. Il£, 20 have studied, as a function of the magnetic field, the behavior of the singlet-triplet splitting 
of the H2-QDM by diagonalizing the two-electron hamiltonian inside the minimal four-dimensional basis 
formed by the products </>L(ri)</>L(r 2 ), </>l(i"i )</>.??, (r 2 ), 4>r(ti)4>l(t2), 0fl(ri)^(r 2 ) of the Is orbitals of 
the separated QD's. This Hubbard- type method (l^] (as well as the refinement employed by Ref. ^(J of 
enlarging the minimal two-electron basis to include the p orbitals of the separated QD's) is an improvement 
over the simple HL method (see Ref. ) , but apparently it is only appropriate for the weak-coupling regime 
at sufficiently large distances and/or interdot barriers. In addition this method fails explicitly (it yields 
a triplet ground state at B — |l9|) for small values of K. Our method is free of such limitations, since 
we employ here an interdot-distance adjustable basis (see section 2) of at least 70 spatial TCO molecular 
orbitals when solving for the sS-UHF ones. Even with consideration of the symmetries, this amounts to 
calculating a large number of two-body Coulomb matrix elements, of the order of 10 6 . 
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UHF orbitals which are self-consistently optimized for any separation d, potential barrier 
height Vj), and magnetic field B. As a result, expression (^) can be characterized as a 
Generalized Valence BondP] (GVB) wave function. Taking into account the normalization 
of the spatial part, we arrive at the following improved wave function for the singlet state 
exhibiting all the symmetries of the original many-body hamiltonian (here, the spatial 
reflection symmetry is automatically restored along with the spin symmetry), 

tf s GV ^(l, 2) = n s V2P V UHF (l, 2) , (10) 

where the normalization constant is given by 

n 2 s = 1/(1 + S UV S VU ) , (11) 

S uv being the overlap integral of the u(r) and v(r) orbitals (which are complex in the 
presence of a magnetic field B), 



S uv = j u*(r)v(r)dr . (12) 
The total energy of the GVB state is given by 

^gvb = n s\^uu ~l~ h-w ~l~ S uv h vu + S vu h uv + J uv + K uv \ , (13) 

where h is the single-particle part (1) of the total hamiltonian (3), and J and K are the 
direct and exchange matrix elements associated with the e — e repulsion e 2 : / 'fir 12. For 
comparison, we give also here the corresponding expression for the total energy of the HF 
"singlet" (i.e., the spin-contaminated determinant with S z = 0), either in the RHF (v = u) 
or sS-UHF case, 

Efjp = h uu + h vv + J uv . (14) 

For the triplet with S z = ±1, the projected wave function coincides with the original 
HF determinant, so that the corresponding energies in all three approximation levels are 
equal, i.e., E t GVB = E^jj F = E\j HF . 



4.3 Comparison of RHF, sS-UHF and GVB results 

In this section, we study in detail the behavior of the interdot coupling in the H 2 -QDM 
in the presence of a perpendicular magnetic field. We present numerical results for three 

5 More precisely our GVB method belongs to a class of Projection Techniques known as Variation before 
Projection, unlike the familiar in chemistry GVB method of Goddard and coworkers [Goddard III, W.A.; 
Dunning, Jr., T.H.; Hunt, W.J.; Hay, P.J. Acc. Chem. Res. 1973, 6, 368], which is a Variation after 
Projection (see Ref. §§). 
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values of the dielectric constant, namely, k = 45 (e — e repulsion much weaker than the case 
of GaAs), k = 25 (e — e repulsion weaker than the GaAs case), and k = 12.9 (case of GaAs). 
In particular we study the evolution of the energy difference, Ae = E s — E t , between the 
singlet and the triplet states as a function of an increasing magnetic field varying from 
B = to B = 9 T. The evolution of the two occupied sS-UHF orbitals of the singlet state is 
illustrated by plotting them at the two end values, B = and B = 9 T. In all three figures 
(i.e., Fig. 2, Fig. 3, and Fig. 4 below), the lower half corresponds to a vanishing interdot 
barrier V& — (unified deformed dot at B — 0), while the upper half corresponds to a finite 
value of Vb, being thus closer to the notion of a molecule proper. The interdot distance is 
chosen to be d = 30 nm in all three cases. 

4.3.1 Case of k = 45 

Fig. 2 displays the evolution of Ae as a function of the magnetic field for k — 45 and for all 
three approximation levels, i.e., the RHF (MO Theory, top solid line), the sS-UHF (dashed 
line), and the GVB (lower solid line). (The same convention for the As(B) curves is followed 
throughout this manuscript.) The case of Vj, — is shown at the bottom panel, while the 
case of Vb = 3.71 meV is displayed in the top panel. The insets display as a function of B 
the overlap (modulus square, l^^l 2 ) of the two orbitals u(r) and v(r) of the singlet state. 
In this calculation, the effective g factor was set equal to zero, g* = 0, so that the gain 
of energy due to the Zeeman effect does not obstruct for large B the convergence of Ae 
towards zero. Due to its smallness relative to hQ e ff = fr(uol + uol/ A) 1 / 2 (where u c = eB/m*c 
is the cyclotron frequency), the actual Zeeman contribution can simply be added to the 
result calculated for g* = 0. 

We observe first that as a function of B, for both the V b — and the V b = 3.71 meV cases 
and for all three levels of approximation, the Ae energy difference starts from a minimum 
negative value (singlet ground state) and progressively increases to zero; after crossing the 
zero value, it remains positive (triplet ground state). However, for large values of B there is 
a sharp contrast in the behavior of the RHF curves compared to the sS-UHF and GVB ones. 
Indeed after crossing the zero axis, the RHF curves incorrectly continue to rise sharply and 
very early they move outside the range of values plotted here (at B — 9 T, the RHF values 
are 0.93 and 1.21 meV for V b = and 3.71 meV, respectively). In contrast, after reaching a 
broad maximum, the positive Ae branches of both the sS-UHF and GVB curves converge 
to zero for sufficiently large values of B. 

The convergence of the singlet and triplet total energies to the same value indicates that 
the H2-QDM dissociates as B attains sufficiently large values. This is also reflected in the 
behavior of the overlaps (see insets) as a function of the magnetic field. In fact, the overlaps 
decrease practically to zero as a function of B, suggesting that the two corresponding 
orbitals u(r) and v(r) of the singelt state tend to become strongly localized on the individual 
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FIGURE 2 



Figure 2: Lateral H2-QDM in the presence of a magnetic field and for k = 45: The left 
column displays the energy difference Ae between the singlet and triplet states according to 
the RHF (MO Theory, upper solid line), the sS-UHF (dashed line), and the GVB approach 
(Projection Method, lower solid line) as a function of the applied magnetic field B. The two 
other columns display the sS-UHF spin- up (f) and spin-down (I) occupied orbitals (modulus 
square) of the singlet state for B = (field-free case, middle column) and B = 9T (right 
column). The top half corresponds to a bare interdot barrier of = 3.71 meV, while 
the bottom half describes the no barrier case V& = 0. For B = 9 T complete dissociation 
has been practically reached. The choice of the remaining parameters is: m* = 0.067m e , 
fojjo = 5 meV, d = 30 nm, and g* — 0. Distances are in nm and the orbital densities in 
1CT 4 nm -2 . Insets: The overlap integral (modulus square) of the two orbitals of the singlet 
state as a function of B. 
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dots. 

The molecular dissociation induced by the magnetic field and the associated electron 
localization is further demonstrated by an examination of the orbital densities (modulus 
square) themselves. These densities are plotted for both the end values of B = (field-free 
case, middle column) and B = 9 T (strong magnetic field, right column). For B = 9 T, it is 
apparent that the plots portray orbitals well localized on the individual dots. In contrast, 
for B = the orbitals are clearly delocalized over the entire QDM. In particular, for B = 
and Vf, = the two orbitals u(r) and v(r), which are different in general, have collapsed to 
the same ls-type distorted orbital associated with the single-particle picture of a unified 
deformed dot. 

One effect of choosing a very weak e — e repulsion (k = 45) is that in the Vj, = case the 
three levels of approximation collapse to the same MO value (no symmetry breaking) for 
magnetic fields below B < 2.9 T, a fact that is also reflected in the orbitals themselves (see 
the B = case in the middle column of the lower half of Fig. 2). However, for V& = 3.71 
meV, even such a weak e — e repulsion does not suffice to inhibit symmetry breaking in the 
field-free case. As a result, for B = and V& = 3.71 meV, the two «(r) and v(r) orbitals, 
although spread out over the entire molecule, are clearly different and the GVB singlet lies 
lower in energy than the corresponding MO value. 

A second observation is that both the sS-UHF and the GVB solutions describe the 
dissociation limit (Ae — > 0) for sufficiently large B rather well. In addition in both the 
sS-UHF and the GVB methods the singlet state at B = remains the ground state for all 
values of the interdot barrier. Between the two singlets, the GVB one is always the lowest, 
and as a result the GVB method presents an improvement over the sS-UHF method both 
at the level of symmetry preservation and the level of energetics. Furthermore, whether a 
singlet or triplet, the GVB always results in a stabilization of the ground state; the improved 
behavior of the GVB over the sS-UHF holds for all values of k. 

We note that the failure of the MO (RHF) approximation to describe the dissociation 
process induced by the magnetic field is similar to its failure to describe dissociation of the 
molecule in the field-free function of the interdot barrier and distance (see Ref. 

2"T|). As discussed in Ref. PH , due to the underlying MO picture, the spin-density func- 



tional calculations of Refs. jT7L [TJ| also fail to describe the molecular dissociation process 
in the field-free case.f] Such spin-density functional calculations are expected to fail in the 

6 Symmetry breaking in coupled QD's within the LSD has been explored by Kolehmainen, J.; Reimann, 
S.M.; Koskinen, M.; Manninen, M. Eur. Phys. J. D 2000, 13, 731. However, unlike the HF case for which a 
fully developed theory for the restoration of symmetries has long been established (see, e.g., Ref. 29 ), the 
breaking of space symmetry within the spin-dependent density functional theory poses a serious dilemma 
[Perdew, J. P.; Savin, A.; Burke, K. Phys. Rev. A 1995, 51, 4531]. This dilemma has not been fully 
resolved todate; several remedies (like Projection, ensembles, etc.) are being proposed, but none of them 
appears to be completely devoid of inconsistencies [Savin, A. In Recent Developments and Applications 
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presence of a magnetic field as well. 



4.3.2 Case of k = 25 

The influence of the magnetic field on the properties of the H 2 -QDM in the case of a stronger 
e — e repulsion (re = 25) is described in Fig. 3. Again we set g* = for the same reasons 
as in the previous case; the meaning of Ae and of the displayed orbitals is the same as in 
Fig. 2. Compared to the k = 45 case, the increase of the e — e repulsion is accompanied by 
an overall strengthening of symmetry breaking and electron localization, since there is no 
range of parameters for which the sS-UHF and GVB solutions collapse into the symmetry- 
adapted MO (RHF) one. Indeed, even for V& = the three levels of approximation provide 
different values for the singlet-triplet difference Ae, with the GVB solution yielding in the 
field- free case the highest stabilization for the singlet. For V& = 3.71 meV, the MO solution 
fails outright even in the field-free case for which it predicts a positive Ae (triplet ground 



state) in violation of the theorem stating |34j] that the ground state of a 2e system at B = 
is always a singlet. 

The overall strengthening of electron localization is also reflected in the behavior of the 
overlap integrals (see insets), since now they start at B = with smaller values compared 
to the corresponding values of the k = 45 case. In addition, the trend towards stronger 
electron localization with smaller k can be further confirmed by an inspection of the orbitals 
of the singlet state displayed in the middle (B = 0) and right (B = 9 T) columns of Fig. 
3. Compared to the corresponding orbitals in Fig. 2, this trend is obvious and we will not 
describe it in detail. It will suffice to stress only that, whatever the degree of initial electron 
localization at B = 0, the applied magnetic field enhances even further this localization, 
achieving again a practically complete dissociation of the artificial molecule at B = 9 T. 



of Modern Density Functional Theory, edited by Seminario, J.M.; Elsevier: Amsterdam, 1996, p. 327]. 
In addition, due to the unphysical self-interaction error, the density-functional theory is more resistant 
against symmetry breaking [see Bauernschmitt, R.; Ahlrichs, R. J. Chem. Phys. 1996, 104, 9047] than the 
sS-UHF, and thus it fails to describe a whole class of broken symmetries involving electron localization, 
e.g., the formation at B = of Wigner molecules in QD's (see footnote 7 in Ref. JlOj] ) , the hole trapping 
at Al impurities in silica [Laegsgaard, J.; Stokbro, K. Phys. Rev. Lett. 2001, 86, 2834], or the interaction 
driven localization-delocalization transition in d- and /- electron systems, like Plutonium [Savrasov, S.Y.; 
Kotliar, G.; Abrahams, E. Nature 2001, 410, 793]. 
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FIGURE 3 



Figure 3: Lateral H2-QDM in the presence of a magnetic field and for k = 25: The left 
column displays the energy difference Ae between the singlet and triplet states according to 
the RHF (MO Theory, upper solid line), the sS-UHF (dashed line), and the GVB approach 
(Projection Method, lower solid line) as a function of the applied magnetic field B. The two 
other columns display the sS-UHF spin- up (f) and spin-down (I) occupied orbitals (modulus 
square) of the singlet state for B = (field-free case, middle column) and B = 9T (right 
column). The top half describes the case of a bare interdot barrier of Vj, = 3.71 meV, while 
the bottom half describes the no barrier case V& = 0. For B = 9 T complete dissociation 
has been practically reached. The choice of the remaining parameters is: m* = 0.067m e , 
fojjo = 5 meV, d = 30 nm, and g* — 0. Distances are in nm and the orbital densities in 
1CT 4 nm -2 . Insets: The overlap integral (modulus square) of the two orbitals of the singlet 
state as a function of B. 
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4.3.3 Case of k = 12.9 (GaAs) 

The dielectric constant for a GaAs heterointerface is k = 12.9, which corresponds to a 
further increase in the e — e repulsion compared to the two previous cases of k — 45 and 
k = 25. This case is presented in Fig. 4. Note that here the Zeeman contribution has 
been included from the beginning by taken g* to be equal to its actual value in the GaAs 
heterointerfaces, i.e., g* = —0.44. As a result, the sS-UHF and GVB curves for Ae converge 
to a straight line representing the Zeeman linear dependence 'jB (with 7 pa 0.026 meV/T 
for S z = +1, short dashed line) instead of vanishing for large B. 

The results presented in Fig. 4 confirm again the trend that electron localization be- 
comes stronger the stronger the interelectron repulsion. As we found previously, the proper 
description of this strong electron localization requires consideration of symmetry breaking 
via the sS-UHF (long dashed curve in the Ae vs. B plots) and the subsequent construction 
of GVB wave functions via Projection Techniques (solid curve in the Ae vs. B plots). The 
MO description is outright wrong, since the corresponding Ae values are positive in the 
whole interval < B < 9 T; in fact they are so large that the whole MO curves lie outside 
the plotted ranges in Fig. 4. For B = 9 T, and in both the barrierless (Vj, = 0) case and 
the case with an interdot barrier of Vj, = 4.94 meV, the artificial molecule has practically 
dissociated. Compared to the previous cases of k — 45 and k = 25, the molecule starts 
with a stronger electron localization already in the field-free case and with increasing B it 
moves much faster towards complete dissociation. 

4.3.4 Overview of common trends in the singlet-triplet energy difference 

According to the GVB calculations, the common trend that can be seen in all cases (inde- 
pendent of the value of k) is that the singlet-triplet energy difference is initially negative 
(singlet ground state) in the field-free case. With increasing magnetic field, this energy dif- 
ference becomes larger (diminishes in absolute value), crosses the value of zero at a certain 
B , and remains positive (triplet ground state) for all B > B , converging from above to 
the straight line representing the Zeeman contribution (or to zero if the Zeeman contri- 
bution is neglected). This means that for sufficiently large values of B the singlet-triplet 
energy difference is given simply by the Zeeman energy and that the molecule has practi- 
cally dissociated. The trend towards dissociation is reached faster the stronger the e — e 
repulsion. 

It is interesting to note again that for k — 45 (weaker e — e repulsion) the sS-UHF and 
GVB solutions collapse to the MO solution for values of the magnetic field smaller than 
B < 2.9 meV. However, for the stronger e — e repulsions (k = 25 and k = 12.9) the sS-UHF 
and GVB solutions remain energetically well below the MO solution. Since the separation 
considered here (d = 30 nm) is a rather moderate one (compared to the value l = 28.50 
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FIGURE 4 



Figure 4: Lateral H2-QDM in the presence of a magnetic field and for k = 12.9 (GaAs): 
The left column displays the energy difference Ae between the singlet and triplet states 
according to the sS-UHF (long dashed line) and the GVB approach (Projection Method, 
solid line) as a function of the applied magnetic field B. The two other columns display the 
sS-UHF spin- up (f) and spin-down (I) occupied orbitals (modulus square) of the singlet 
state for B = (field- free case, middle column) and B = 9 T (right column). The top 
half describes the case of a bare interdot barrier of V& = 4.94 meV, while the bottom half 
describes the no barrier case V& = 0. For B = 9 T complete dissociation has been practically 
reached. The choice of the remaining parameters is: m* = 0.067m e , TujOq = 5 meV, d — 30 
nm, and g* = —0.44 (GaAs). Distances are in nm and the orbital densities in 10~ 4 nm~ 2 . 
Insets: The overlap integral (modulus square) of the two orbitals of the singlet state as a 
function of B. 
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nm at B = for the extent of the Is lowest orbital of an individual dot with Hujq = 5 meV), 
we conclude that there is a large range of materials parameters, interdot distances, and 
magnetic-field values for which the QDM's are weakly coupled and cannot be described by 
the MO theory; a similar conclusion was reached in Ref. [|ST| for the field-free case. 



5 Conclusions 

We have shown that, even in the presence of an applied magnetic field, the sS-UHF method, 
in conjunction with the companion step of the restoration of symmetries when such sym- 
metries are broken, is able to describe the full range of couplings in a QDM, from the 
strong-coupling regime exhibiting delocalized molecular orbitals to the weak-coupling one 
associated with Heitler-London-type combinations of atomic orbitals. 

The breaking of space symmetry within the sS-UHF method is necessary in order to 
properly describe the weak-coupling and dissociation regimes of QDM's. The breaking 
of the space symmetry produces optimized atomic-like orbitals localized on each individ- 
ual dot. Further improvement is achieved with the help of Projection Techniques which 
restore the broken symmetries and yield multideterminantal many-body wave functions. 
The method of the restoration of symmetry was explicitly illustrated for the case of the 
H2-QDM in the presence of a magnetic field. It led to the introduction of a Generalized 
Valence Bond many-body wave function as the appropriate vehicle for the description of 
the weak-coupling and dissociation regimes of artificial molecules. 

Additionally, we showed that the RHF, whose orbitals preserve the space symmetries 
and are delocalized over the whole molecule, is naturally associated with the molecular 
orbital theory. In a generalization of the field-free case of natural [36—38] and artificial [ET 



molecules, it was found that the RHF fails to describe the weak-coupling and dissociation 
regimes of QDM's in the presence of an applied magnetic field as well. 



This research is supported by a grant from the U.S. Department of Energy (Grant No. 
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